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ABSTRACT 

CGRaBS J0809+5341, a high redshift blazar at z = 2.144, underwent a giant opti¬ 
cal outburst on 2014 April 19 when it brightened by ~5 mag and reached an unfiltered 
apparent magnitude of 15.7 mag. This implies an absolute magnitude of —30.5 mag, 
making it one of the brightest quasars in the Universe. This optical flaring triggered 
us to carry out observations during the decaying part of the flare covering a wide en¬ 
ergy range using the Nuclear Spectroscopic Telescope Array, Swift, and ground based 
optical facilities. For the first time, the source is detected in 7 -rays by the Large 
Area Telescope onboard the Fermi Gamma-Ray Space Telescope. A high optical po¬ 
larization of ~10% is also observed. Using the Sloan Digital Sky Survey spectrum, 
accretion disk luminosity and black hole mass are estimated as 1.5 x 10 45 erg s _1 and 
10 8 ' 4 M 0 respectively. Using a single zone leptonic emission model, we reproduce the 
spectral energy distribution of the source during the flaring activity. This analysis sug¬ 
gests that the emission region is probably located outside the broad line region, and the 
jet becomes radiatively efficient. We also show that the overall properties of CGRaBS 
J0809+5341 seems not to be in agreement with the general properties observed in high 
redshift blazars up to now. 

Subject headings: galaxies: active — gamma rays: galaxies — quasars: individual 
(CGRaBS J0809+5341) — galaxies: jets 
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1. Introduction 

Blazars are a subclass of active galactic nuclei (AGN) with powerful relativistic jets aligned 
close to the line of sight to the observer (Urry & Padovani 1995). They emit over the entire 
electromagnetic spectrum, predominantly by non-thermal emission processes. Because of small 
inclination angle, the emission from their jet is Doppler boosted. Blazars are classified as flat 
spectrum radio quasars (FSRQs) and BL Lac objects based on the rest frame equivalent width 
(EW) of their broad optical emission lines with FSRQs having EW > 5 A (Stocke et al. 1991; 
Stickel et al. 1991). However, Ghisellini et al. (2011a) (see also Ghisellini et al. 2009) have 
recently proposed a new classification scheme based on the luminosity of broad line region (BLR) 
measured in units of Eddington luminosity with FSRQs having Lblr/LecM > 5 x 10 -4 . Both 
classes share many common properties, such as flat radio spectra ( a r < 0.5; S u oc v~ a ) at 
GHz frequencies, rapid flux and polarization variations (Wagner & Witzel 1995; Andruchow et al. 
2005) and exhibit superluminal patterns at radio wavelengths (Jorstad et al. 2005). 

The broadband spectral energy distribution (SED) of blazars consists of two distinct peaks. 
The low energy peak lies between the regimes of infrared (IR) and X-ray while the high energy 
one lies in the MeV—TeV range. The low energy component is known to result from synchrotron 
emission whereas the origin of the high energy component is still a matter of debate. In the 
leptonic emission model, the high energy peak in the SED is explained by the inverse-Compton 
(IC) scattering of synchrotron photons from the jet (synchrotron self Compton or SSC, Konigl 
1981; Marscher & Gear 1985; Ghisellini & Maraschi 1989). Alternatively, the seed photons for 
IC scattering can be external to the jet (external Compton or EC, Begelman & Sikora 1987; 
Melia & Konigl 1989; Dermer et al. 1992). The high energy component can also be explained as 
a result of hadronic processes (see e.g., Miicke et al. 2003; Bottcher et al. 2013). FSRQs and BL 
Lac objects are found to follow the so-called blazar sequence (Fossati et al. 1998; Ghisellini et al. 
1998). However, recently Giommi et al. (2012) have proposed that the existence of such sequence 
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could be due to selection effects. 

CGRaBS J0809+5341 (hereafter J0809+5341) is a high redshift FSRQ (z = 2.144 ± 0.002, 
Paris et al. 2014) which was overlooked for a long time due to its faintness and/or prolonged 
quiescence. It was predicted as a candidate 7 -ray emitter by Healey et al. (2008), but was not 
detected by any earlier 7 -ray surveys. It is radio bright (F 8A ghz = 183.4 mJy, Healey et al. 2007), 
having a flat radio spectrum and exhibits compact core morphology in the Faint Images of the 
Radio Sky at Twenty centimeters (FIRST) observations. 

In this work, we present a detailed multi-wavelength study of J0809+5341 which recently 
flared in the optical band. By analyzing both archival and simultaneous observations, we present 
a consistent picture of its emission processes in the framework of leptonic radiation models 
of blazars in its low and high activity states. In Section 2, we describe the multi-wavelength 
campaign carried out to study the optical outburst of this source. The details of the data reduction 
procedure are reported in Section 3 and the results are presented in Section 4. We discuss our 
findings in Section 5 and provide conclusions in Section 6 . Throughout the work, we adopt a 
ACDM cosmology with the Hubble constant H 0 = 71 km s _1 Mpc~\ Q m = 0.27, and Ha = 0.73. 


2. Optical Outburst And Multi-Wavelength Campaign 

J0809+5341 underwent a giant optical outburst on 2014 April 19 (Balanutsa et al. 2014), 
when its unfiltered apparent magnitude reached 15.7 mag (see Figure 1). This corresponds to 
an absolute magnitude of —30.5 mag, thereby making it one of the brightest quasars in the 
Universe. Compared to the archival Sloan Digital Sky Survey (SDSS) observations, the source 
was ~ 5 mag brighter on 2014 April 19. In order to study this extraordinary event, we organized a 
multi-frequency campaign using both space and ground based telescopes. The hard X-ray mission 
Nuclear Spectroscopic Telescope Array ( NuSTAR , Harrison et al. 2013) observed the target on 
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2014 May 8 and was supplemented by a simultaneous Swift target of opportunity (ToO) and 
optical polarimetric observations from the Telescopio Nazionale Galileo 1 (TNG). Additionally, 
the source was observed on 2014 April 26, 27, and 28 by other Swift ToO observations. We 
observed the source in Bessel B, V, R, and I filters on 2014 April 27 and in Bessel U, B, V, R, 
and I bands on 2014 May 1 from the Himalayan Chandra Telescope 1 (HCT). The entire campaign 
was complemented by continuous monitoring from the Large Area Telescope (LAT) onboard the 
Fermi Gamma-ray Space Telescope (hereafter Fenni-LAT). 


3. Multiwavelength observations and Data Reduction 
3.1. Fermi- Large Area Telescope Observations 

The Fermi- LAT data used in this work were collected over the first 72 months of 
Fermi operation, from 2008 August 5 (MJD 54,683) to 2014 August 4 (MJD 56,873). We 
follow the standard data analysis procedures as mentioned in the Fermi-LAI documentation * 2 3 . 
Events belonging to the SOURCE class and in the energy range 0.1 —300 GeV are used. A filter 
“DATA_QUAL>0” && “LAT_CONFIG==l” is used to select good time intervals and a cut of 
100° is applied on the zenith angle to avoid contamination from the Earth limb 7 -rays. Recently 
released galactic diffuse emission component glljem_v05_revl.fit and an isotropic component 
iso_source_v05_revl.txt are considered as background models 4 . The normalization parameter of 
the background models are kept free during the fitting. The binned likelihood method included 
in the pylikelihood library of Science Tools (v9r33p0) along with the post-launch 

'httpi/Avww.tng.iac.es 

2 http: //www. iiap .res. in/centers/iao 

3 http://fermi. gsfc.nasa.gov/ssc/data/analysis/documentation/ 

4 http://fermi. gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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instrument response functions P7REP_SOURCE_V15 are used in the analysis. 

The significance of the 7 -ray signal is evaluated by means of a maximum likelihood test 
statistic TS = 2A log(£) where C represents the likelihood function, between models with and 
without a point source at the position of the source of interest. Sources lying within 10° region 
of interest (ROI) centered at the position of J0809+5341 and defined in the second Fermi -LAT 
catalog (2FGL; Nolan et al. 2012) are included in the model file. The spectral parameters of all 
the nearby sources are taken from the 2FGL catalog and are allowed to vary (except the scaling 
factor) during the likelihood fitting. In addition to that, we also include the sources lying between 
10° to 15° from the center of the ROI and keep their parameters fixed to the 2FGF values. 
Further, Fermi-LAY has detected many sources after the release of the 2FGF catalog and hence 
these sources are not included in it 5 . If lying close to the source of interest, these unmodeled 
sources could affect the results of the analysis. The presence of these sources, if any, is tested 
by generating the residual TS map of the ROI. We find two new sources (Figure 2), whose 
positions are optimized using the tool gtfindsrc and obtained as, R.A., Dec. = 122°.136, 
49°.737 (J2000) and 132°.506, 51°. 157 (J2000) respectively. The possible counterparts of these 
7 -ray sources should be the FSRQ OJ 508 and a 7 -ray emitting narrow line Seyfert 1 galaxy 
SBS 0846+513 respectively (see The Fermi-FAT Collaboration 2015). We verify that no other 
significant sources are left in the data to be modeled, by generating another residual TS map 
including these two new sources in the model file. We do not find any new sources with TS > 25 
(see Figure 2). The two new sources are then modeled with a power law model and included in 
the analysis. 

A first round of the likelihood fitting is performed over the entire 72 months of the FAT data 
and all the sources with TS < 25 are removed from the model. This updated model is then used 
for further time series and spectral analysis. To generate the spectra and the light curves we allow 


5 http://www.cbpf.br/icrc2013/papers/icrc2013-1153 .pdf 
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the normalization parameter of all the sources within the ROI to vary, whereas the photon indices 
are fixed to the values obtained from the average analysis over appropriate time intervals. Light 
curves are generated by adopting the unbinned likelihood method as it is expected to encounter 
fewer events over shorter time intervals 6 . The source is considered to be detected if TS > 9 which 
corresponds to ~ 3a detection (Mattox et al. 1996). For 1 < TS < 9, we calculate 2a upper limit. 
This is done by varying the flux of the source till TS reaches a value of 4 (Abdo et al. 2010). 

We do not calculate upper limits if TS < 1. Primarily governed by uncertainty in the effective 
area, the measured fluxes have energy dependent systematic uncertainties of around 10% below 
100 MeV, decreasing linearly in log(E) to 5% in the range between 316 MeV and 10 GeV and 
increasing linearly in log(E) up to 15% at 1 TeV 7 . Errors associated with the LAT data analysis 
are the la statistical uncertainties, unless otherwise specified. 


3.2. NuSTAR Observations 

The source was observed with NuSTAR for a cleaned exposure time of 28 ks. The NuSTAR 
data are reduced and filtered for background flares using the NuSTAR Data Analysis Software 
(NUSTARDAS) version 1.4.1, and response files are generated using CALDB version 20140814 
Spectra are extracted for the two focal plane modules (FPMA and FPMB) using the nuproducts 
tool. The source spectrum is extracted from a 30" circular region centered on the peak emission, 
and the background spectrum is extracted from a 70" circular region on the same chip, free of 
contaminating sources. The spectra are source dominated until > 20 keV, and they are binned to 
have a minimum of 20 counts per bin. 

6 http://fermi. gsfc.nasa.gov/ssc/data/analysis/scitools/likelihood _tutorial.html 

7 http://fermi. gsfc.nasa.gov/ssc/data/analysis/LATxaveats.html 
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3.3. Swift Observations 

The source is below the sensitivity limit of the Burst Alert Telescope (BAT; Barthelmy et al. 
2005). However, it is significantly detected by the pointed observations from the X-Ray Telescope 
(XRT; Burrows et al. 2005) and the Ultraviolet Optical Telescope (UVOT; Roming et al. 2005). 
This is the first time that J0809+5341 is detected in the X-ray band. 

The XRT data are processed with the XRTDAS software package (v.3.0.0) available within 
the HEASOFT package (6.16). Event files are cleaned and calibrated using standard filtering 
criteria with xrtpipeline (v. 0.13.1 ) and the calibration database that was updated on 
2014 July 30. Standard grade selections of 0-12 in the photon counting mode are used. Cleaned 
event files are then summed using the task XSELECT. To extract the source spectrum from 
the summed event files, a circular region of 20 pixel (~47") centered at the source position, is 
chosen, while background is extracted from a nearby circular region of 50 pixel radius. All the 
exposure maps are combined with XIMAGE and used to generate ancillary response files using 
the task xrtmkarf . The source spectrum is binned to have at least 1 count per bin. Spectral 
fitting is done using XSPEC (Arnaud 1996). An absorbed power law (Ah = 3.75 x 10 20 cm^ 2 ; 
Kalberla et al. 2005) is used for fitting. Due to the faintness of the source, we use C-statistics 
(Cash 1979) within XSPEC. The uncertainties are calculated at 90% confidence level. 

UVOT has observed the source in 4 filters namely V, B, U, and UVW1. All the observations 
are integrated using the task uvotimsum and analyzed with uvot source. Source region is 
selected as a circle of 5" radius centered at the source position. Background is chosen from 
a nearby source free circular region of 1' radius. The observed magnitudes are corrected for 
reddening using the galactic extinction of Schlafly & Finkbeiner (2011) and converted to flux 
units using the zero point magnitudes of Breeveld et al. (201 1). The optical-UV flux of the high 
redshift blazars (z > 2) could be significantly absorbed by neutral Hydrogen in the intervening 
Lyman a absorption systems. We use the attenuation calculated by Ghisellini et al. (2010b) to 
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correct for this effect. 


3.4. Ground Based Optical Observations 

3.4.1. Polarimetry 

Linear polarimetry was carried out on 2014 May 8 at the TNG. Observations were performed 
with the PAOLO polarimeter 8 with the Sloan r filter for about an hour. Data reduction as well 
as aperture photometry are performed using custom-made software tools 9 . Instrumental Stokes 
parameters, polarization degree, and position angle are then corrected for the instrumental 
polarization by means of a suitable number of polarized and un-polarized standard stars as 
described in detail in Covino et al. (2014). Flux calibration is derived by the observation of the 
SDSS isolated and non-saturated stars in the polarimeter field of view (see next section). No flux 
variability, total or polarized, is detected during the observations. 


3.4.2. Photometry 

As mentioned earlier, J0809+5341 was observed using HCT on two epochs, namely 2014 
April 27 in B, V, R, and I filters and on 2014 May 1 in U,B,V, R, and I filters. The details of 
the instrument can be found in Paliya et al. (2014). Standard procedures in IRAF are used to do 
the pre-processing of the images (bias subtraction, flat-fielding, and cosmic ray removal). After 
pre-processing, instrumental magnitudes of the target as well as stars in the field are obtained via 
PSF photometry using DAOPHOT (Stetson 2011) available in MIDAS 10 . We could not convert 

8 http://www. tng.iac.es/instruments/lrs/paolo.html 

9 https://pypi.python.org/pypi/SRPAstro.FITS 

1() Munich Image Data Analysis System 
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these instrumental magnitudes to standard magnitudes directly as Landolt standard star fields 
were not observed during the observations. Also, there are no stars in the field of the source with 
UBVRI photometry available. Therefore, to convert the measured magnitudes to the standard 
system the following procedure is adopted. We obtain ugriz magnitudes for 3 comparison stars 
in the field from the SDSS 11 (SDSS J080938.95+533813.8, SDSS J080956.48+534341.3 and 
SDSS J081006.16+534248.0). For these three stars the UBVRI magnitudes are derived using the 


following transformation equations (Jordi et al. 2006) 

U = B + (0.52 ± 0.06) *(u-g) + (0.53 ± 0.09) * (g - r) - (0.82 ± 0.04) (1) 

B = g+ (0.313 ± 0.003) * (g - r) + (0.219 ± 0.002) (2) 

V = 1+ (0.675 ± 0.002) *(g-i) + (0.364 ± 0.002) (3) 

R = I+ (0.930 ± 0.005) * (r - i) + (0.259 ± 0.002) (4) 

I = i- (0.386 ± 0.004) *(i- z )~ (0.397 ± 0.001) (5) 


Once the UBVRI magnitudes of these three stars are obtained, the derived instrumental 
magnitudes of J0809+5341 are converted to UBVRI magnitudes using differential photometry. 
The standard magnitudes are then dereddened and converted to flux using the zero points of 
Bessell et al. (1998). The flux are corrected for absorption by intervening Lyman a absorption 
systems as done by Ghisellini et al. (2010b). 

Recently, Carrasco et al. (2014) have reported the NIR observations of J0809+5341 from 
the 2.1 m telescope of the Guillermo Haro Observatory operated by the National Institute for 
Astrophysics, Optics and Electronics (Mexico). We take their reported magnitudes in J, H, and Ks 
bands, correct for reddening and convert to flux units using the zero points of Bessell et al. (1998). 

n http://skyserver.sdss3.org/drl0 
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4. Results 


4.1. Black Hole Mass And Accretion Disk Luminosity 


J0809+5341 has a low S/N spectrum from SDSS-DR10 (Paris et al. 2014). As this is 
the only spectrum available for this source, we use it to estimate the black hole mass (Mbh) 
and the accretion disk luminosity (La The spectrum is brought to the rest-frame and then 
dereddend using E(B—V) of 0.038 taken from the NED 12 . The fitting to the SDSS spectrum is 
based on chi-square minimization using the MPFIT package (Markwardt 2009). A power law 
continuum is fit to the spectrum using the line free regions [1445, 1465] A and [1700,1710] A, 
i.e. on either side of the C IV line. This continuum is then subtracted from the spectrum. Fe 
emission is not subtracted from the spectrum as it is known to be weak for C IV line (Shen et al. 
2011 ). The modified spectrum, between the wavelength range [1500, 1600] A, is then fit with 
a single Gaussian function. A narrow component C IV is not considered in our fitting as virial 
black hole mass estimates using C IV are based on the FWHM of the entire C IV line profile 
(Vestergaard & Peterson 2006). The fit to the C IV line is shown in Figure 3. From the single 
Gaussian fit we find the line flux and the a of the C IV line as (46.46 ± 0.76) x 10" 1 ' erg cm" 2 
s" 1 A " 1 and 1338.3 ± 139 km s" 1 respectively. Correcting the observed a for the resolution of 
the instrument, the FWHM of C IV is estimated as 3145 ± 327 km s" 1 . Using the region between 
1345 — 1350 A, we find the mean continuum flux at 1350 A to be (4.96 ± 0.26) x 10" 17 erg cm" 2 
s" 1 A" 1 . This translates to a continuum luminosity (XL\) of (2.43 ± 0.14) x 10 45 erg s" 1 . The 
C IV line luminosity is obtained as Lq i v = (1.69 ± 0.03) x 10 43 erg s" 1 . Using the measured 
FWHM and the continuum luminosity, we estimate the black hole mass using the following 
equation (Shen et al. 201 1) 


log 


/ -^BH,vir\ 

V m q ) 


a + b log 


XL} 


10 44 ergs 



+ 2 log 


/FWHM\ 
\ kms" 1 / 


( 6 ) 


12 http://ned.ipac.caltech.edu/ 
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The adopted values of a and b are 0.66 and 0.53 respectively, taken from Shen et al. (2011). This 
results in a black hole mass of log ( 777 ^) = 8.39 ± 0.21. It is important to note that for jet 
dominated sources like blazars, single epoch black hole mass estimates may not be relied upon, 
particularly in the flaring states as (i) the BLR will be ionized by non-thermal continuum emission 
and (ii) the BLR cannot be considered to be in equilibrium during episodes of strong flaring 
activity (Leon-Tavares et al. 2013). However, since the available SDSS spectrum corresponds 
to the faint state of the source, the black hole mass estimated here seems robust. Following 
Celotti et al. (1997) (see also Francis et al. 1991), we calculate the total BLR luminosity (L B lr) 
using the flux of the C IV line. We find L BL r = 1.5 x 10 44 erg s _1 . Assuming 10% of the accretion 
disk luminosity is reprocessed by the BLR, the accretion disk luminosity is 1.5 x 10 45 erg s _1 . 


4.2. Average Gamma-Ray Properties 

J0809+5341 is not present in the 2FGL catalog, indicating that it was not detected with TS 
> 25 in first two years of Fermi operation. Indeed, the LAT data analysis for this period results 
in TS = 0.5, explaining the absence of the source in the 2FGL catalog. However, analysis of 
the third through fifth year of the LAT data gives TS ~ 185.5 (~ 13a, Mattox et al. 1996), thus 
confirming that J0809+5341 is a high redshift 7 -ray emitting FSRQ. Therefore, for the first time, 
we report the detection of J0809+5341 in the 7 -ray band 13 . Moreover, the source is found to 
be in a relatively bright state during the sixth year of Fermi operation with TS = 191.3 and the 
derived 7 -ray flux and photon index are (2.78 ± 0.44) x 10 -8 ph cnr 2 s" 1 and 2.15 ± 0.08 
respectively. The details of the results of the average analysis of Fermi -LAT data, covering various 
time intervals, are given in Table 1. 

13 This source is now included in the recently released 3FGL catalog as 3FGL J0809.5+5342 
(The Fermi-LAT Collaboration 2015). 
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We perform 7 -ray point source localization over the photons extracted during the third 
through fifth year of Fermi operation (the time period when the source was detected significantly) 
using the tool gt f indsrc. The localized spatial coordinates are RA = 122°.446, Dec = 53°.673 
(J2000), at an angular separation of 0°.02 from the radio position of J0809+5341 (RA = 122°.424, 
Dec = 53°.690, J2000), with a 95% error circle radius of 0°.04. This implies a close spatial 
association of the 7 -ray source with the radio counterpart. 

In order to search for the energy of the highest energy photon, we use the tool gt srcprob 
and event class CLEAN. The highest energy is found to be 13.63 GeV, detected on 2013 December 
31, at 0°.03 far from the source position, and with 98.4% probability to be associated with the 
7 -ray source. 


4.3. Gamma-ray Temporal Variability 

The 7 -ray light curve of J0809+5341, covering the first 72 months of Fermi operation, is 
shown in Figure 4. Fermi-LAI data are binned monthly and we also show daily scaled 7 -ray flux 
covering the period of high optical activity (MJD 56,748—56,786 or 2014 April 1 to 2014 May 
9), by blue circles in the inset. It is clear from Figure 4 that the source was not detected at all in 
the first two years of Fermi operation. It was detected by the FAT sporadically during the third 
through fifth year at a low flux level. However, it becomes relatively active only during the last 
12 months when it was continuously detected by the FAT. Visual inspection of the daily binned 
light curve shown in the inset of Figure 4 hints for the presence of two peaks approximately at the 
same flux level (one at around MJD 56,760 and other at MJD 56,775). Moreover, we calculate 
the daily binned 7 -ray flux (F 7 ) for the flaring period and find a maximum F 7 of (2.56 ± 0.96) x 
10 - ' ph cm -2 s -1 on 2014 April 14 (MJD 56,761), which coincides with the first report of the 
high optical activity from the source (Shumkov et al. 2014). 
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4.4. Optical-UV Observations 

The results of multi epoch Swift -UVOT observations are given in Table 2. It is evident 
from this table that the source has shown flux variations on both a daily and weekly timescales. 
Also, multi-band apparent magnitudes, as observed from HCT, are presented in Table 3. For a 
comparison with the archival data, we also give the SDSS magnitudes converted to UBVRI filters 
using the transformations given by Jordi et al. (2006). These data are not corrected for Galactic 
reddening. It is clear from Table 3 that the source has brightened significantly in all the filters 
compared to the archival observations. 

The polarimetric observations from TNG on 2014 May 8 show the detection of high optical 
polarization from J0809+5341. The observed polarization is found to be as high as 9.8 ± 0.5%, 
which is typically seen in blazars. The corresponding observed polarization angle is 98 ± 1 deg. 


4.5. Spectral Analysis 

We test the presence of curvature in the overall y-ray spectrum of the source by fitting a 
LogParabola model. It is defined as dN/dE oc (E/E 0 )~ a ~^ l ° 9< ^ E I E °\ where E 0 is an arbitrary 
reference energy fixed at 300 MeV, a is the photon index at E 0 and /3 is the curvature index 
which defines the curvature around the peak. The test statistic of curvature is then evaluated as 
T Scu rv e = 2 (log £ (LogParabola) — log £(power-law)). We set the threshold TS CUIve > 16 to 
test the presence of a significant curvature as done by Nolan et al. (2012). There are hints of the 
presence of curvature as TS curve ~ 10 (~ 3cr; Table 1), though a strong claim can not be made as 
TS CUIve is below the threshold of 16 set in the 2FGL catalog (Nolan et al. 2012). 

The joint XKY/NuSTAR spectrum is well fit with a simple power law (C-statistic/d.o.f.= 
40/35), modified by Galactic absorption, and we find a best-fit photon index of 1.4 ± 0.1. This 
fit is shown in Figure 5, along with the residuals to the model. We allow for a difference in flux 
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calibration between NuSTAR and XRT spectra by including a constant (CONST in XSPEC) fixed 
at 1 for two NuSTAR spectra (calibration difference between FPMA and FPMB are on the order 
of 1%, and are not detectable in spectra with few counts, see e.g., Marinucci et al. 2014) and 
free to vary for the XRT. This constant offset is consistent with 1 (1.2 ± 0.3) and we find the 
same value of photon index whether or not we allow for this constant multiplicative factor. The 
large error is caused by the low count rate, and low sensitivity of the XRT compared to other 
X-ray telescopes. We calculate the goodness of the fit using Monte-Carlo simulations with the 
GOODNESS command in XSPEC, and find that 47% of the simulated spectra based on the model 
have a lower y 2 . 

Fitting the absorbed power law model to the combined XRT spectra of the first three 
observations on 2014 April 26, 27 and, 28 (for a total exposure of 4.8 ksec yielding 48 counts) 
resulted in a photon index of 1.2 ± 0.4 while that derived for the fourth observation on 2014 
May 8 (for a total exposure of 4.6 ksec yielding 40 counts) is 1.9 ± 0.5. This suggests a possible 
softening of the X-ray spectrum over a course of ~10 days. However, considering the large errors 
in the photon indices, due to the faintness of the source, a firm conclusion cannot be reached. 

The B—R color of the source, from HCT observations, are found to be 0.90 ± 0.04 and 0.94 
± 0.05 for the epochs of 2014 April 27 and 2014 May 1 respectively and thus there is no color 
change between these two epochs. At the epoch of the SDSS observation, we find B—R = 0.98 ± 
0.06. Thus, we do not see any optical color variation between the epochs 2003 and 2014, though 
the source has varied significantly in flux. 


4.6. Spectral Energy Distributions 


We generate the SED of the source during two different activity states. A high activity 
state covering the period of the recent optical outburst (2014 April 1 to 2014 May 9), and a low 


- 16- 


activity state for which the LAT data covering the third through fifth year of Fermi operation 
(see Figure 4) and other non-simultaneous archival observations 14 are used. The derived flux 
values are given in Table 4. Assumption of the archival observations as low activity state can 
be justified by the fact that during the recent optical flare, the optical magnitude was brighter by 
~ 5 mag compared to the archival SDSS measurements (Balanutsa et al. 2014). However, since 
the SDSS observations were taken well before the launch of Fermi satellite, the results based on 
such non-simultaneous data can be questioned. The source was observed by Mobile Astronomical 
System of the TElescope-Robots (MASTER; Lipunov et al. 2010) on 2011 March 25, i.e. between 
the third to fifth year of Fermi operation (see Figure 1). The obtained upper limit in the unfiltered 
magnitude was 19.4 mag. This hints that during the third through fifth year of Fermi operation, 
the optical flux level of the source was possibly similar to that observed by the SDSS in 2003. 
Moreover, as can be seen in Figure 4, the source become active only very recently, therefore the 
third through fifth year of Fermi-LAI observations can be adopted as low activity state. 

We use a simple one zone leptonic emission model, similar to the one used by 
Ghisellini & Tavecchio (2009, hereafter GT09) and Finke et al. (2008) (see also Dermer et al. 
2009), to interpret the SEDs of the source. The emission region is assumed to be spherical, filled 
with relativistic electrons, located at a distance Zdiss, from the central black hole of mass 
and moving relativistically with bulk Lorentz factor F. A scaling factor of T oc Z 0.5 

is assumed 

(where Z is the distance from the black hole) in the inner parts of the jet where it is anticipated to 
be accelerating and parabolic in shape (GT09, Vlahakis & Konigl 2004). After the acceleration 
phase, the jet becomes conical with semi aperture angle 0.1 rad and a constant bulk Lorentz factor 
r 

1 max* 

The electron energy distribution is assumed to follow a smoothly joining broken power law of 
the form N'(Y) = Ajj ^yyy^y/y ^ ,where p and q are the particle indices before and after the 

14 http ://www. asdc. asi.it/bzcat 
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break energy ( 7 ^) respectively (primed quantities are measured in the comoving frame). The size 
of the emission region is adopted by considering it to cover the entire jet cross-section. Thermal 
emission from the accretion disk is evaluated assuming a standard optically thick, geometrically 
thin disk (Shakura & Sunyaev 1973) with inner and outer radii R in = 3i?s c h, R 0 u t = 500f?s c h 
respectively (GT09), where Rs c h is the Schwarzschild radius. Locally, the accretion disk spectrum 
is explained by a multi-temperature blackbody (e.g. Frank et al. 2002). Above and below the 
accretion disk, the presence of the X-ray corona is also considered which reprocesses 30% of 
the accretion disk luminosity. The inner and outer radii of the corona are assumed to be 3f?sch 
and 30i?sch respectively. The spectrum emitted by the corona is adopted as cut-off power law: 
L cor (e) oc e“" cor exp(—e/e c ) (GT09), where e is the dimensionless photon energy (= 77 ^). The 
cut-off energy is assumed to be 150 keV and we adopt a flat power law i.e. cc COI = 1- The BLR 
is assumed as a spherical shell located at a distance X B lr = lO 17 ^^ cm, where L d ,45 is the 
accretion disk luminosity in units of 10 45 erg s ” 1 (GT09). It reprocesses 10% of the accretion 
disk luminosity and its spectrum is assumed to be a blackbody peaking at rest-frame Lyman-a 
frequency (Tavecchio & Ghisellini 2008). The dusty torus, for simplicity, is assumed to be a 
thin spherical shell, located at a distance of Z m = 10 18 cm and reprocess 50% of the 
accretion disk emission. Its emission profile can be explained by a simple blackbody peaking 
at temperature T IR . In the presence of randomly oriented magnetic field, electrons radiate via 
synchrotron and IC scattering mechanisms. The synchrotron and synchrotron self-Compton 
emissions, in the observer’s frame, are calculated using the formulations of Finke et al. (2008). 
Besides these mechanisms, the particle also loose energy through IC scattering of external photon 
field from the accretion disk (EC-disk), BLR (EC-BLR), and the dusty torus (EC-torus) (GT09; 
Dermer et al. 2009; Dermer & Menon 2009; Cao & Wang 2013). Finally, the kinetic power of the 
jet is calculated assuming both protons and electrons to have equal number densities (Celotti et al. 
1997). Protons are assumed to be cold and contribute only to the inertia of the jet. For the SED 
modeling, we adopt the black hole mass as 10 8A M & and the accretion disk luminosity as 1.5 x 
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10 45 erg s _1 , as calculated in Section 4.1. Since the SED in the low activity state is generated 
using non-simultaneous data, we perform the modeling only on the flaring state SED where the 
source was monitored contemporaneously over a wide energy range. In Figure 6, we show the 
model spectrum due to different emission mechanisms along with the observed fluxes and the 
corresponding parameters are given in Table 5. The variation of the radiation energy densities, 
measured in the comoving frame, are also shown in the bottom panel of Figure 6. To model the 
SEDs, we start with a plausible set of parameters which are then constrained within a range that 
better represent the data. 


5. Discussion 

The non-detection of J0809+5341 by earlier high energy missions indicates that either the 
source is intrinsically faint or remains in quiescence for a long time. The recent optical outburst 
along with the increased brightness across the electromagnetic spectrum made it possible to detect 
the source for the first time in X-rays. This also led to the first detection of the source in the y-ray 
band as predicted earlier by Healey et al. (2008). 

In Figure 1, we show the MASTER image of J0809+5341 taken on 2011 March 25 and 2014 
April 19. As can be seen clearly, the source was not detected in 2011 and only an upper limit in 
the unfiltered magnitude was obtained. However, it becomes extremely bright in 2014. The y-ray 
light curve also shows that J0809+5341 became active only in the sixth year of Fermi operation 
when enhanced y-ray emission is observed during the optical outburst. Comparing the results of 
the y-ray analysis done from the third through fifth year of Fermi operations with that obtained 
during the sixth year reveals that (i) the flux increased in the sixth year, and (ii) the photon index 
does not change within the errors. 


During the period of high activity, the maximum one day binned y-ray photon flux is found 
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to be (2.56 ± 0.96) x 10^' ph cm -2 s -1 which corresponds to an isotropic 7 -ray luminosity (L 7 ) 
of 9.3 x 10 48 erg s” 1 . This, in turn, corresponds to a luminosity measured in the proper frame of 
the jet as L 7 em ~ L 7 /2T 2 ~ 1.2 x 10 46 erg s _1 considering a bulk Lorentz factor T = 20 obtained 
from the SED modeling (Table 5). This is a good fraction of the kinetic jet power (~ 53%; P-j.kin 
= 2.2 x 10 46 erg s _1 ) indicating that the jet becomes radiatively efficient and a significant amount 
of the kinetic jet power gets converted to the radiative power. 

The SEDs of J0809+5341 reveal that the optical-UV spectrum is steep which we interpret 
as synchrotron emission. Alternatively, one can associate this spectrum to the accretion disk 
emission. However, it is unlikely that the large variation seen at optical-UV energies during 
low and high activity states can be a result of the perturbations in the accretion disk emission. 
Moreover, the analysis of the SDSS spectrum also suggests a relatively less luminous disk 
(Section 4.1). In addition, the source has shown high optical polarization during the recent flare, 
thus supporting the synchrotron origin of the optical-UV spectrum. The X-ray spectrum of the 
source, obtained from the 5wz//-XRT and NuSTAR observations during the flare, is typical of 
powerful blazars and can be well reproduced by the SSC process. 

It can be seen in Figure 6 that the shape of the 7 -ray spectrum is relatively flat compared 
to the optical-UV spectrum. A steep optical-UV spectrum suggests that the spectral shape of 
the electron energy distribution is soft. If electrons of similar energies are contributing to the 
emission at 7 -rays through IC process, then one would expect a steep 7 -ray spectrum similarly 
to the optical-UV. However, we observed a hard 7 -ray spectrum. Hence, we assume that the 
harder 7 -ray spectrum is the result of the superposition of different IC emission processes. One 
possibility is that the SSC is contributing in the hard X-ray to soft 7 -rays and an EC component 
explaining the remaining spectra, however, this is not possible since the steep falling optical 
spectrum suggests that the synchrotron emission peaks at lower energies which in turn causes 
negligible SSC emission at high energies (e.g., Sahayanathan & Godambe 2012). This indicates 
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that, though significant, SSC cannot explain the observed 7 -ray spectrum. Alternatively, a flat 
7 -ray spectrum can be a result of interplay between various EC mechanisms. We find that the 
inclusion of EC-BLR and EC-torus emission can reproduce the 7 -ray spectrum satisfactorily. 
The relative contribution of these target photons measured in the comoving frame depends on the 
location of the emission region from the central engine. Then, under the assumption that the 7 -ray 
spectrum is the superposition of EC-BLR and EC-torus components, we can derive the location 
of the emission region by reproducing the observed 7 -ray spectral shape. In the bottom panel of 
Figure 6 , we plot the energy densities of various components in the comoving frame, as a function 
of Z diss . Modeling the flaring SED indicates that the emission region is located at a distance 
from the central black hole where both the BLR and IR-torus energy densities are contributing 
almost equally to produce the observed 7 -ray spectrum. We further constrain the parameters of 
the present model by considering near-equipartition between relativistic particles and magnetic 
field. The resulting model spectrum along with the observed fluxes are shown in Figure 6 and the 
corresponding parameters are given in Table 5. 

Comparison of the flaring state SED with that representing the low activity state indicates 
few interesting features. The increase of the optical flux appears to be higher than that of the 7 -ray 
flux. This is supported by the fact that, at the time of the flare, the source was one of the brightest 
quasars in the optical band whereas the rise of the 7 -ray flux was relatively modest. However, it 
should be noted that we do not have simultaneous 7 -ray observations at the time of low activity 
state as recorded by the SDSS optical monitoring and thus a strong claim by comparing the 
non-simultaneous observations cannot be made. Further, comparison of the accretion disk flux 
with that of the archival optical observations hints that even during quiescence the optical-UV 
emission is synchrotron dominated. 

The SEDs and associated modeling parameters of J0809+5341 are quite different from those 
obtained for other high redshift blazars. Generally, the optical-UV part of the SED of high redshift 
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blazars is dominated by extremely luminous accretion disk radiation, the peak of IC emission lies 
in the hard X-ray regime resulting in a steep 7 -ray spectrum, and they also are known to host 
more than a billion solar mass black hole at their centers (e.g., Ghisellini et al. 2011b, 2013). 

In contrast, J0809+5341 hosts a relatively less luminous accretion disk and less massive central 
black hole. The optical-UV spectrum is dominated by synchrotron radiation, and the IC peak lies 
at GeV range. Thus, the overall observed properties of J0809+5341 indicates that this source is, 
in many ways, different from other high redshift blazars but show similarities with its low redshift 
counterparts (Ghisellini et al. 2010a). 


6. Conclusions 

In this work, we present a detailed multi-frequency study of the high redshift blazar 
J0809+5341. The main findings of the work are summarized below. 

1. Predicted as a candidate 7 -ray emitter by Healey et al. (2008), J0809+5341 is now detected 
in the 7 -ray band by Fermi-LAI, as confirmed by the 3FGL catalog. 

2. The black hole mass and accretion disk luminosity, calculated by analyzing the archival 
SDSS spectrum, are found to be 1O 8 - 1 2 3 4 5 M 0 and 1.5 x 10 45 erg s _1 respectively. 

3. A significant fraction of the kinetic jet power gets converted to the radiative power during 
the flare, or in other words, the jet becomes radiatively efficient. 

4. The flaring state optical-UV spectrum can be successfully modeled by the synchrotron 
emission. The choice of the synchrotron mechanism over the accretion disk is primarily 
influenced by the observation of high optical polarization during the flare and by the flare 
itself. 

5. The optical-UV spectrum is found to be steep and the 7 -ray spectrum is relatively flat. The 
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flatness of the 7 -ray spectrum can be explained by locating the emission region outside the 
BLR where both the BLR and torus energy densities play a major role in describing the 
observed 7 -ray spectrum. 

6 . Many of the observed properties of J0809+5341 are at odds with that generally observed in 
other high redshift blazars (Ghisellini et al. 201 lb, 2013). 
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Table 1: Details of the power law model fits to the averaged 7 -ray data for various time periods. The 
quoted flux values are in units of 10 ~ 8 ph cm -2 s -1 whereas L 1 is the 7 -ray luminosity. The last 
column quotes the significance of the curvature present in the spectrum by means of a LogParabola 
model fitting. 


Time Period (MJD) 

To. 1-300 GeV 

-^0.1—300 GeV 

TS 

log L 7 

TQ 
u curve 

54683-55412 (2008 Aug 05 to 2010 Aug 4) 

- 

- 

0.5 

- 

- 

55412-56508 (2010 Aug 04 to 2013 Aug 4) 

2.26 ± 0.08 

1.61 ±0.26 

185.5 

47.72 

8.70 

56508-56873 (2013 Aug 04 to 2014 Aug 4) 

2.15 ±0.08 

2.78 ± 0.44 

191.3 

48.00 

9.64 

54683-56873 (2008 Aug 05 to 2014 Aug 4) 

2.28 ± 0.06 

1.40 ±0.19 

267.6 

47.65 

9.60 


Table 2: Swift -UVOT observations of J0809+5341. The quoted numbers are in magnitudes. 


UVOT-filters 

2014 April 26 

2014 April 27 

2014 April 28 

2014 May 8 

UVW1 

- 

19.17 ±0.12 

- 

18.62 ± 0.07 

U 

17.39 ±0.08 

- 

17.93 ± 0.09 

17.28 ± 0.04 

B 

18.30 ±0.14 

- 

18.69 ±0.11 

18.00 ±0.05 

V 

17.58 ±0.05 

- 

17.97 ±0.13 

17.42 ± 0.07 
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Table 3: Apparent magnitudes of J0809+5341, as observed from HCT. For comparison, the ugriz 
magnitudes from SDSS that are converted to UBVRI magnitudes are also given. 


Filters 2014 April 27 2014 May 1 2003 November 20 (SDSS) 


u 


- 


18.15 

± 

0.07 

21.03 

± 

0.12 

B 

18.14 

± 

0.04 

18.48 

± 

0.04 

21.19 

± 

0.04 

V 

17.67 

± 

0.02 

17.99 

± 

0.03 

20.61 

± 

0.04 

R 

17.17 

± 

0.01 

17.47 

± 

0.03 

20.21 

± 

0.04 

I 

16.56 

± 

0.02 

16.82 

± 

0.04 

19.69 

± 

0.05 
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Table 4: Summary of the SED analysis 


Fermi-L AT 

Activity state a 

Period b 

F{). 1 -300 GeV C 

P d 

1 0.1-300 GeV 

TS e 

^pred 

Q 

55412-56508 

1.61 ±0.26 

2.26 ± 0.08 

185.47 

897.63 

F 

56748-56786 

5.87 ± 1.19 

2.14 ±0.08 

50.48 

110.94 

NuSTAR 

Activity state 

Exp. g 

F h 

J- 3-79 keV 

^3-79 keV 1 

Normal ization' 

Stat. k 

F 

28 

1.30 ±0.30 

1 S o+0.65 
1-50—0.46 

3.30 ± 0.40 

12/14 

Swift-X RT 

Activity state 

Exp. g 

ro.3-10 keV 1 3 4 5 6 7 

Eb.3-10 keV™ 

Normalization' 

Stat k 

F 

9.45 

1 4<;+0.30 
1 - hj -0.29 

4 S 1+1-32 

5 09 +1,29 

52.14/74 

Swift- UVOT 

Activity state 

V n 

B n 

U n 

UVWl n 


F 

1.99 ±0.09 

1.56 ±0.09 

1.64 ± 0.06 

0.62 ± 0.02 



1 Different activity states selected for the modeling; Q: low activity and F: high activity state. 

2 Time period considered for the SED modeling, in MJD. 

3 Integrated 7 -ray flux in 0.1—300 GeV energy range in units of 10 ~ 8 ph cm ' 2 s" 1 . 

4 Photon index calculated from 7 -ray analysis. 

5 Significance of detection using likelihood analysis. 

6 Number of predicted photons during the time period under consideration. 

7 Net exposure in kiloseconds. 

8 Photon index of the power law model. 

9 Observed flux in units of 10 -12 erg cm -2 s _1 , in 3—79 keV energy band. 

^Normalization at 1 keV in 10 ~ 5 ph cm -2 s _1 keV” 1 . 

Statistical parameters: C-stat./dof. 

12 Photon index of the absorbed power law model. 

13 Unabsorbed flux in units of 10 -13 erg cm -2 s _1 , in 0.3 — 10 keV energy band. 

14 Average flux in Swift V, B, U, and UVW1 bands, in units of 10 “ 12 13 14 erg cm -2 s _1 . 
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Table 5: Summary of the parameters used/derived from the modeling of the SED. 


Parameter 

Symbol 

High Activity 

Slope of particle spectral index before break energy 

P 

1.6 

Slope of particle spectral index after break energy 

<7 

5.3 

Magnetic field in Gauss 

B 

2.0 

Particle energy density in erg cm~ 3 

u ' 

0.05 

Bulk Lorentz factor 

r 

20 

Minimum Lorentz factor 

^min 

1 

Break Lorentz factor 

ib 

1714 

Maximum Lorentz factor 

*1max 

le5 

Size of the BLR in the units of f?sch 

Zblr 

1570 

Distance of the emission region from the black hole in parsec ( Rs c h ) 

■^diss 

0.048(1980) 

Dusty torus temperature in Kelvin 

Eir 

900 

Viewing angle in degrees 

0 

3 

Jet power in electrons in log scale 

Pe 

44.62 

Jet power in magnetic field in log scale 

Pb 

45.11 

Radiative jet power in log scale 

Pr 

46.19 

Jet power in protons in log scale 

Pp 

46.34 
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Fig. 1.— Mobile Astronomical System of the TElescope-Robots (MASTER) image of 
J0809+5341 (field of view of 5" x 5") taken on 2011 March 25 (left) and on 2014 April 19 (right). 
An upper limit of 19.4 mag is obtained during 2011 observations, while there is ~ 5 mag brighten¬ 
ing (with respect to the SDSS observations) during the 2014 optical outburst. North is up and East 
is left. 
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Fig. 2.— Left: Residual TS map of the 0.1 —300 GeV events during the time period covered in this 
work, centered on the coordinates of J0809+5341. Two new sources in the ROI are noticed. Right: 
Residual TS map of the same region after modeling both the new sources. 
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Rest-frame Wavelength (A) 


Fig. 3.— Top: A fit to the continuum and C IV line in the SDSS spectrum of J0809+5341. Bottom: 
The residual between the fit and the observed spectrum. 
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Fig. 4.— Flux history of J0809+5341 in the 7 -ray band, covering the first 72 months of Fermi- 
LAT operation. Data are in units of 10 -8 ph cm -2 s^ 1 . 95% upper limits are shown by downward 
arrows. Q represents the low activity state, while F denotes the high activity state selected for the 
SED modeling. Inset: Daily binned 7 -ray light curve covering the period of high activity. The flux 
units are same as that of the main panel. 
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Fig. 5.— Top: Joint XRT (from 0.3-7 keV) and NuSTAR (from 3-60 keY) spectrum, fit with an 
absorbed power-law. Bottom: data to model ratio for this fit. Data are rebinned for clarity, and the 
two NuSTAR spectra are grouped in Xspec but fit separately. 
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Fig. 6.— Top left: Low activity state SED of J0809+5341. Gray bow-tie plot represents third 
through fifth year of Fermi operation. Black downward arrow is the first two years 95% flux upper 
limit, calculated by assuming the photon index obtained from the analysis of the third through fifth 
year LAT data. Top right: Modeled high activity state SED. Black dotted line represents thermal 
contributions from the torus, accretion disk, and X-ray corona. Pink thin solid line and green long 
dashed line are synchrotron and SSC radiation. Grey dashed, red dash-dot, and black dash-dot-dot 
lines represent EC-disk, EC-BLR, and EC-torus components respectively. Blue thick solid line 
refers to the sum of all the radiation components. Bottom panel: Variation of the energy densities 
measured in the comoving frame, as a function of the distance from the black hole, in units of f? Sch . 
Vertical line denotes the location of the emission region. 
















